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Summary 

This  report  summarizes  the  results  of  a  one  year  basic  research  effort  under  the  Air  Force  Office 
of  Sponsored  Research  contract  No.  F49620-92-J-0450 

As  summarized  in  our  proposal  and  letter  to  Dr.  Spencer  Wu,  the  stated  goal  of  this  research  has 
been 

“...During  the  first  year  of  our  research  we  will  derive  the  analytical  model 
of  the  aeroservoelastic  system  to  account  for  the  coupled  fluid  flow/control/ 
structure  field  equations.  Our  methodology  to  accomplish  this  research  goal 
above  will  be  to  integrate 

(1)  the  infinite  dimensional,  game  theoretic  control  theory  for  distributed 
parameter  systems  developed  by  [Kurdila]  and, 

(2)  a  model  of  fluid  flow  /  structure  interaction  as  developed  by  [Strganac]. 

It  is  anticipated  that  we  will  additionally  further  identify  issues  to  complement  the 
needs  of  the  research  community.” 

Our  research  accomplishments  have  far  exceeded  our  stated  goals  for  the  first  year.  During  this 
period, 

(1)  we  have  extended  the  infinite  dimensional,  game  theoretic  control 
strategy  to  account  for  structured  uncertainty  arising  in  the  control  influence 
operator  associated  with  piezoceramic  actuators, 

(2)  we  have  extended  the  infinite  dimensional,  game  theoretic  control 

strategy  to  account  for  structured  uncertainty  arising  in  constitutive  law  evolution, 

(3)  we  have  derived  a  coupled  fluid/structure/control  analysis  that  synthesizes 
the  distributed  control  approach  derived  in  (i)  and  (ii)  above  with  a  flutter  induced 
phenomenological,  microcrack  model  for  a  panel. 

Our  first  year  of  research  has  not  only  accomplished  the  stated  goals  of  the  proposal,  but  has  also 
enabled  the  research  team 

(1)  to  develop  an  experimental  wind-tunnel  facility  designed  expressly  for 
the  study  of  the  onset  of  flutter,  nonlinear  aeroelasticity  and  verification  of  the 
control  theories  derived  in  this  research; 

(2)  to  derive  multiresolution-based  analysis  methods  for  the  investigation 
of  the  onset  of  flutter  in  our  experiment^  facilities; 


During  pilot  studies  conducted  under  AFOSR  contract  F49620-92-J-0450  it  has  been 
demonstrated  that  wavelet-Galerkin  methods  are  directly  applicable  to  a  wide  class  of  fluid-struc¬ 
ture-control  interaction  problems  of  relevance  to  the  Air  Force.  For  examples  figures  (3  111) 
through  (3.1.1-5)  are  attracted  extracted  from  [KKKS]  and  illustrate  that  wavelet-Galerkin  meth¬ 
ods  provide  excellent  representations  of  sharp  edge  effects  commonly  introduced  in  modelline 
piezoceramie  actuators.  Figure  (3.1. 1-1)  depicts  the  control  influence  function  represented  via  a 
Daubechies  order  3  wavelet,  and  the  structured  uncertainty  function  employed  to  represent  edee 
delamination.  Figure  (3. 1.1 -2)  illustrates  disturbance  attenuation  in  a  closed-loop  response 
denved  using  an  infinite  dimensionel,  game-theoretic  control  design  employing  a  wavelet  basis. 
Figures  (3. 1.1 -3)  and  (3. 1.1-4)  illustrate  the  improved  performance  of  the  game  theoretic  control 
m  comp^son  to  linear  quadratic  regulation.  Both  control  simulations  are  conducted  using  the 
wavelet-Galerkin  method.  Figure  (3.1.1-5)  illustrates  the  profound  difference  between  finite  and 
mfimte  dimensional  stability  margins  calculated  when  employing  the  wavelet-Galerkin  method. 

other  studies  supported  by  AFOSR  contract 
F49620-92-J-{H50  It  has  become  clear  that  a  central  issue  in  employing  multiresolution  analysis 
in  theoretical  mechanics  and  dynamics  is  the  derivation  of  general  and  efficient  techniques  for 

Ilac  approaches  have  been  derived  under  AFOSR  con¬ 

tract  F49620-92-J-0450,  and  will  be  presented  at  the  35th  Structures,  Structural  Dynamics  and 
Matenals  Conference  m  April.  1994.  These  three  methods  are  based  upon  the  use  of 

(i)  boundary  integral  equations, 

(ii)  numerical  boundary  measures,  and 

(iii)  wavelet-based  finite  elements 

to  derive  w^el^-based  representations  of  fluid/structure/control  interaction  problems.  For  exam- 
p  e,  figure  (3.1-6)  illustrates  the  use  of  boundary  integral  equations  to  represent  fluid  flow  around 
a  two  dimensional  airfoil  using  raultiresolulion  analysis  methods.  These  results  are  extracted 
irom  IKulJ  and  have  been  carried  out  under  AFOSR  contract  F49620-92-J-()450. 

.  of  representing  coupled  fluid-structure  interaction  problems  is  pre¬ 

sented  in  [KKKS]  and  [KKW],  wherein  numerical  boundary  measures  are  utilized  to  enable  the 
so  ution  of  fluid-structure  interaction  problems  on  general  domains.  Figures  (3. 1.1 -7)  and  (3.1.1- 
8)  enable  the  visuali^tion  of  a  simple  airfoil  and  its  associated  numerical  boundary  measure. 
Figures p.  1.1 -9)  and  (3.1.1-12)  illustrate  the  central  conclusion  of  [KKW]. 

It  is  possible  to  generate  numerically  stable  boundary  measures  that  retain  the 

nonnegative  character  of  the  underlying  boundary  value  problem.” 

Specifically,  this  conclusion  is  apparent  from  the  results  of  “method  2”  in  figures  (3  1  1-91 
through  0.1.1-12).  The  interested  reviewer  i  referred  to  [KKW]  for  a  more  detailed  discussion 
Figures  (31.1-13)  through  (3.1.1-14)  demonstrate  that  it  is  also  possible  to  combine  nonlinear 
dymmical  system  theory  with  computational  mechanics  to  solve  multiresolution  problems  on 
^bitr^  domains.  In  figure  (3.1.1-13),  an  iterated  function  system  is  constructed  such  that  the 
boundary  of  interest  is  its  domain  of  attraction.  Figure  (3.1.1-14)  illustrates  the  domain  embed¬ 
ding  solution  obtained  employing  an  associated  ergodic  theory  described  in  [Ku2].  As  noted  ear- 


displacement  (nondimensionalized) 


Figure  (3.1. 1-1)  Control  Influence  and  Disturbance 


Figure  (3. 1.1-2)  Closed  Loop  Response  to  Initial  Condition 


Figure  (3.1. 1-3)  Disturbance  Attenuation,  Game 
Theoretic  Compared  to  Linear  Quadratic  Regualtion 
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Figure  (3. 1.1-4)  Disturbance  Rejection,  Game 
Theoretic  Compared  to  Linear  (Quadratic  Regualtion 


Figure  (3. 1.1 -5)  Robust  Stability  Discrepancy,  Finite  Dimensional 
Models  Compared  to  Infinite  Dimensional  Limit 


Figure  (c)  Fast  Wavelet  Transform  of 
Solution,  Level  1  Resolution 


Figure  (e)  Multiresolution  Representation, 
3  Levels,  Overwritten,  Daubechies  Order  3 


Figure  (d)  Fast  Wavelet  Transform  of  Solution, 
Level  3  Resolution 


Figure  (0  Wavelet  Galerkin  Solution  of 
Boundary  Integral  Equation,  potential  function 


Figure  (3. 1.1-6)  Wavelet-Galerkin  Solution  of  Boundary  Integral 
Equations  for  Flow  Past  An  Airfoil,  Multiresolution  Representation 


Minimum  Eigenvalue 


Figure  (3. 1.1 -9)  Negative  Eigenvalue  Magnitude 
as  a  Function  of  Penalty  Parameter 


Figure  (3.1.1-10)  Number  of  Negative  Eigenvalues 
as  a  Function  of  Penalty  Parameter 


Figure  (3.1.1-11)  Number  of  Iterations  for  Convergence 
as  a  Function  of  Penalty  Parameter 
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Figure  (3.1.1-12)  Smoothness  of  Domain  Embedding,  Comparison  of  Numerical 
Boundary  Measures,  Method  1  (right)  and  Method  2  (left) 
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Closely  Couple  Airfoils,  (64x64)  Pi  ^  (3.1.1-I8)  Closely  Couple  Airfoils,  (128x128) 

Potential  Function  Contours  Potential  Function  Contours 


lier,  the  length  details  of  this  theory  can  be  found  in  [Ku2]. 

Finally,  the  work  in  [KKP]  and  depicted  in  figures  (3.1.1-15)  through  (3.1.1-18)  demon¬ 
strate  how  the  results  of  [LRT]  and  [DM]  can  be  extended  to  derive  wavelet-based  finite  elements. 
These  specialized  techniques  provide  an  effective  means  for  representing  more  complex  geome¬ 
try,  such  as  the  study  of  flow/structure  interaction  effects  in  closely  coupled  airfoils  depicted  in 
figures  (3.1.1-15)  through  (3.1.1-18). 

QlDevelonment  of  Anahtkal  Flutter  Anahsk  and  Experimentnl  FncilitiP^ 

Aeroelastic  instabilities,  which  may  be  catastrophic,  occur  as  a  result  of  coupling  between 
aerodynamic  forces,  structural  forces,  and  inertial  forces.  Therefore,  the  aeroelastic  design  of  air¬ 
craft  must  address  the  structure  and  the  aerodynamic  forces  within  the  flight  environment.  The 
proposed  program  of  research  will  initially  use  a  structural  experiment  limited  to  two  degrees  of 
freedom  shown  schematically  in  figure  (3.2.2- 1).  The  actual  experimental  setup  is  depicted  in  fig¬ 
ure  (3.2.2-3)  and  has  been  incorporated  in  the  3’  x  4’  subsonic  wind  tunnel  at  Texas  A&M. 

The  structure  has  been  designed  to  implement  a  continuous  nonlinear  response  which  rep¬ 
resents  an  extension  of  recent  research  in  piecewise  linear  models  which  account  for  a  deadband 
in  linear  response.  The  equations  of  motion  for  the  nonlinear  aeroelastic  model  representing  the 
wing  section  limited  to  the  pitch  and  plunge  degrees  of  freedom  can  be  written  as: 

I'a  H-  mrycos  (a  +  (]))  +  (a  +  boc’)  -i-  C^a.  =  M 
my  +  mracos  («  +  ([))-  mrasin  («  +  (())+  (y  +  ay^)  +  C^y  =  -Ncos  (a  +  (])) 

where 

N=N[y, y,  a,  a,  /;([),  AR, planform] 

M  =  M[y,  y,  a,  d,  AR,  planform] 

and  where  the  overdots  represent  time  derivatives,  y  and  a  are  the  plunge  and  pitch  coordinates, 
Cy  and  Cpj  are  the  damping  coefficients  for  the  support,  Ky  and  are  the  spring  constants,  a  and 

b  are  of  nonlinear  spring  constants,  (j)  is  the  static  angle  of  attack,  m  is  the  mass,  /  is  the  mass 

moment  of  inertia  about  the  elastic  axis,  r  is  the  distance  between  the  elastic  axis  and  the  center  of 
the  mass,  N  is  the  aerodynamic  normal  force,  and  M  is  the  aerodynamic  pitching  moment. 

The  linear  form  of  these  equations  are  found  in  several  texts  on  the  subject  of  aeroelastic- 
ity.  The  equations  are  coupled  by  the  motion  and  the  parameter,  r.  The  equations  appear  to  be 
uncoupled  for  the  case  when  r  =  0;  yet,  it  is  important  to  note  that  the  aerodynamics  of  the  above 
equations  can  take  many  forms  but  even  the  most  elementary  models  of  unsteady  aerodynamic 
l^havior  provide  coupling  between  the  two  degrees  of  freedom.  Nonlinear  models  for  the  transla¬ 
tional  and  torsional  springs  have  been  employed  in  these  equations  as  well  as  the  experimental 
setup. 


Simulation  results  for  a  two  degree-of-freedom  model  of  the  actual  structure  which  use  a 
linear  model  of  the  structure  are  presented  in  Figure  3.2.2-2.  The  results  illustrate  classical 
aeroelastic  behavior.  Tj  and  T2  represent  periods  unique  to  the  individual  degrees  of  freedom. 
However,  in  this  case  we  note  a  transition  to  a  period  of  7j.  The  frequencies  have  coalesced  and 
the  motion  in  the  rotational  degree  of  freedom  (alpha)  is  increasing  -  both  indications  of  an 


aeroelastic  instability.  The  complexity  of  the  response  actually  observed  in  the  wind  tunnel  for 
the  experimental  setup  is  evident  from  figures  (3.2.2-4)  and  (3.2.2-5),  which  represent  the  tran¬ 
sient  response  of  the  nonlinear  system  in  pitch  and  plunge,  respectively.  Obviously,  the  response 
of  the  experimental  system,  that  is  subject  to  nonlinearities  in  aerodynamics  and  structural 
dynamics,  differs  profoundly  from  that  predicted  by  linear  theory  in  figure  (3.2.2-2).  Just  as 
importantly,  it  is  extremely  difficult  to  locate  or  predict  when  these  forms  of  pathological  response 
will  occur.  This  observation  reflects  an  important  feature  of  several  of  the  current  methods  of 
study  of  nonlinear  dynamics.  That  is,  many  researchers  concern  themselves  with  the  observation 
and  identification  of  pathological  response  in  a  given  system,  be  it  a  physical  system  or  a  mathe¬ 
matical  model.  The  research  herein  characterizes  the  transition  to  instabilities.  This  latter 
approach  is  of  much  greater  value  to  the  designer.  In  comparison  to  the  integral  wavelet  trans¬ 
forms  of  the  last  section,  figures  (3.2.2-6)  and  (3.2.2-7)  illustrate  the  time-frequency  representa¬ 
tions  that  closely  resemble  the  bifurcation  diagrams  commonly  used  for  the  study  of  instability  in 
nonlinear  dynamical  system  theory.  This  proposal  will  investigate  the  use  of  the  “coalescence  of 
frequencies  with  scale”  evident  in  the  wavelet  transforms  depicted  in  figures  (3.2.2-6)  and  (3.2.2- 
7)  as  a  means  of  characterizing  transition  to  instability. 
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Abstract 

Wavelet  analysis  is  regarded  as  an  extremely  promis¬ 
ing  tool  for  approximate  solution  of  multi-field  prob¬ 
lems,  such  as  those  arising  in  aeroelasticity  and 
fluid  structure  interaction,  due  to  its  inherent  multi¬ 
resolution/multi-scale  nature.  However,  wavelet  analy¬ 
sis  has  been  conducted  primarily  within  the  fields  of 
signal  and  image  processing  due  to  the  difficulty  in 
defining  wavelet  bases  that  satisfy  specified  boundary 
conditions.This  paper  employs  an  embedded  dommn 
technique  to  ameliorate  the  difficulty  associated  with 
deriving  a  wavelet  basis  for  a  specific  multi-field  ini¬ 
tial/boundary  value  problem.  Instead  of  constructing 
an  explicit  wavelet  basis  over  the  domain  of  interest, 
boundary  conditions  are  enforced  using  a  penalty  for¬ 
mulation  that  requires  the  calculation  of  a  numerical 
boundary  measure.  This  paper  presents  strategies  for 
the  rapid  calculation  of  numerical  boundary  measures 
employed  in  wavelet-Galerkin  approximations  of  prob¬ 
lems  in  aeroelastic  transient  response  and  control.  In 
addition,  the  impact  of  new  wavelet  quadrature  trunca¬ 
tion  error  bounds  is  discussed  in  the  context  of  aeroe¬ 
lastic  simulation  and  control. 


(I)  Introduction 

In  1988,  Daubechies  introduced  a  class  of  compactly 
supported  wavelet  functions  which  give  a  new  decompo¬ 
sition  of  Ir^(R),  and  by  tensor  products,  L^(R”).  Since 
these  functions  combine  orthogonality  with  localization 
and  scaling  properties,  there  have  been  a  number  of 
applications  in  engineering,  especially  in  image  and  sig¬ 
nal  processing  But  recently,  the  application  of 

wavelets  in  the  approximate  solution  for  partial  differ¬ 
ential  equation  has  become  an  area  of  active  research. 


For  the  most  part,  this  research  has  focussed  on  two 
distinct  problems; 

(i)  representation  of  arbitrary  boundary  conditions  in  a 
wavelet  basis,  and 

(ii)  derivation  of  wavelet  bases  that  are  tailored  to  the 
differential  operator  under  consideration. 

Representative  work  dealing  with  these  topics  can 
be  found  in  Both  of  these  problems  are  indeed 
formidable,  and  are  guaranteed  to  be  the  source  of  re¬ 
search  over  the  next  few  years.  One  difficulty  associated 
with  the  representation  of  arbitrary  boundary  condi¬ 
tions  is  that  some  wavelet  bases  cannot  be  expressed  in 
closed  form.  Combined  with  the  fact  that  the  wavelet 
bases  are  often  not  interpolatory,  prescribing  basis  sys¬ 
tems  of  desired  smoothness  and  assuming  specified  val¬ 
ues  at  the  boundary  can  be  difficult.  With  respect  to 
the  second  class  of  problems,  some  progress  has  been 
made  for  specific  combinations  of  differential  operators 
and  boundary  conditions.  Within  these  tailored  basis 
systems,  the  differential  operator  can  be  expressed  in  a 
diagonal  form  and  results  in  simplified  computational 
schemes.  As  one  might  expect,  however,  the  class  of  dif¬ 
ferential  operators  and  boundary  conditions  for  which 
this  may  be  accomplished  is  rather  limited  and  progress 
in  research  is  made  slowly. 

This  paper  circumvents  the  difficulties  associated 
with  adapting  a  wavelet  basis  to  a  given  differential  op>- 
erator  or  domain  by  employing  a  penalty  implementa¬ 
tion  of  the  governing  differential  equations.  The  penalty 
formulation  embodies  domsdn  embedding  methods  that 
have  a  well  established  theoretical  foundation  in  com¬ 
putational  mechanics  As  will  be  demonstrated  in 
this  paper,  the  critical  task  in  this  formulation  is  the 
rapid  and  accurate  calculation  of  specific  boundary  in¬ 
tegrals  associated  with  the  boundary  conditions.  These 
integrals  may  be  approximated  using  numerical  bound¬ 
ary  measures  introduced  by  or  by  employing  alterna- 


tlve  approximations  derived  in  the  paper.  Tiie  numer¬ 
ical  comparison  of  the  wavelet  galerkin  formulations  in 
aeroelasticity  is  based  upon  the  panel  flutter  transient 
Tesponse  and  control  problem  studied,  for  example,  in 
^  and 


(II)  The  Fluid-Structure-Control  Model 


Aeroelasticity  is  defined  as  the  phenomenon  result¬ 
ing  from  the  interaction  between  aerodynamic,  inertial, 
and  elastic  forces.  While  the  class  of  problems  in  fluid- 
structure  interaction  to  which  the  theory  presented  in 
this  paper  is  vast,  the  panel  flutter  problem  has  been 
selected  to  demonstrate  the  underlying  principles  in 
wavelet  galerkin  approximation  methods.  The  panel 
flutter  problem  has  been  studied  by  many  distinguished 
researchers  over  the  years  /  0  and  has  recently  re¬ 
ceived  renewed  interest  due  to  advances  in  actuation 
using  active  materials  1.  Dowell  ^2  briefly  describes 
the  change  in  the  aeroelastic  stability  arising  from  the 
degradation  of  isotropic  panels.  In  this  research,  ^2  sug¬ 
gests  that  nonlinear  flutter  analysis  permits  the  predic¬ 
tion  of  the  fatigue  life  of  the  associated  panel  structure. 
This  work  shows  that  the  conventional  fatigue  data  may 
be  used  for  estimating  the  life  of  the  panel  at  flutter  con¬ 
ditions  if  stress  levels  and  frequency  are  known.  Mei  et, 
al.  ^3  studies  the  fatigue  life  of  beams  with  thermal 
effects.  Frequencies  and  stresses  for  limit  cycle  motions 
of  these  structures  subject  to  various  thermal  loadings 
are  determined.  The  critical  dynamic  pressure  and  fa¬ 
tigue  life  history  are  determined.  These  results  provide 
a  dynamic  pressure  at  which  infinite  fatigue  life  is  gu- 
ranteed.  Chen  et.  al.  ^0  present  flutter  studies  of  thin 
cracked  panels.  The  structure  is  assumed  to  behave  in  a 
linear  manner  and  only  single  cracks  are  considered.  A 
liybrid-displacement  finite  element  is  used  in  order  to 
facilitate  modeling  the  singularity  near  the  crack  tip. 
The  influence  of  crack  lengths  and  flow  directions  on 
the  critical  dynamics  pressure  is  presented. 

As  noted  earlier,  several  relatively  new  open  ques¬ 
tions  have  arisen  regarding  the  simulation  and  control 
of  panel  flutter  using  actuation  devices  comprised  of  ac¬ 
tive  materials.  Scott  presents  an  optimization  strat¬ 
egy  based  on  linear  quadratic  control  design  to  suppress 
flutter  using  adaptive  piezoelectric  materials.  While 
presents  an  interesting  approach  to  the  use  of  adaptive 
materials  to  reduce  stress  concentrations  in  the  vicinity 
of  a  crack  tip,  few  authors  have  considered  the  effect 
of  damage,  evolution,  or  loss  of  actuation  authority  in 
control  strategies  employing  adaptive  materials.  The 
current  paper  extends  the  recent  work  in  ^3  on  robust 
control  of  distributed  parameter  systems  by  considering 
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where 

Hi{x)=H{x-ai)-H{x-Pi)  (2) 

and  where  p  is  the  material  density,  h  is  the  laminae 
thickness,  El  is  the  bending  modulus,  ci  is  the  viscous 
damping  coefficient,  ^*1  is  the  stiffness  constant  of  the 
beam  foundation,  and  gi  is  the  piezoelectric  gain.  Each 
gi  depends  upon  the  specific  laminate  geometry.  The 
interested  reader  should  see  M  or  ^5  for  details.  The 
remainder  of  this  paper  will  be  concerned  with  the 
wavelet-galerkin  approximation  of  the  above  equation 
for  obtaining  transient  response  and  control  design. 


(Ill)  Multiresolution  Analysis 


All  derivations  in  this  paper  employ  Daubechies 
wavelets,  which  are  introduced  in  ^6.  For  each  positive 
iV,  there  is  a  family  of  compactly  supported  wavelets  of 
dimension  iNT  —  1.  Let  </>  be  the  Daubechies  scaling  func¬ 
tion  of  order  N{n4>  in  ^6),  and  let  {/loj  *  ‘  ‘  be 

the  corresponding  coefficients  for  the  2-scale  difference 
equation 

2N-1 

^(x)  =  ^  hk<l>i2x  -  k)  (3) 

k=0 


where  the  coefficients  {/iq, /ii, •  •  • ,  h2N^i}  satisfy 


2N-1 


s/2  hk  =  2 

k=0 

(4) 

2N-1 

7^  hkhk-\^2l  =  5o,z,  /  €  Z 

k=0 

(5) 

where  is  the  Kronecker  delta  function.  By  choosing 
:=  {—l)^'^^h2N-i-k  (6) 

One  can  find  compactly  supported  wavelet  functions 
which  are  in  L^(R)  and  which  are  defined  using  a 
fundamental  wavelet 


2;v-i 


fc=0 


(7) 


The  functions  <)>  and  ip  both  have  support  in  the  interval 
[0,  2N  —  1].  The  wavelet  system  associated  with  h  is 
defined  by  the  relationship 

;  (f>{(x)  :=  2^^^(j>{Vx  -  k) 

:=  2^^^'ip{2^x  —  k) 


(8) 


’wlicre  €  Z,  7  E  Z'^(tlie  nouuegative  integers).  Then 
the  above  functions  have  following  properties. 


supp  V'i(x))  = 


k  +  2N  -I 
2i  ’  2-’ 
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Y^<l>i{x)  =  2\  xeR  (11) 

kez 

[  x^x[;{x)dx  —  0,  771  =  O,-- ,  -  1  (12) 

JR 

Moreover  where,  for  example,  a(2)  ^ 

0.55,  a(3)  1.088,  a(4)  «  1.618. 

Daubechies  wavelets  induce  a  specific  multiresolu- 
tion  analysis  of  L^(R).  Let,  for  j  E  Z, 

Vj  :=  closure(span{(?!>{(x) :  k  E  Z})  (13) 

and  let  Wj  be  the  orthogonal  complement  of  Vj  in  Vj^i^ 
then  for  fixed  j  E  Z,  {(/>{.}  is  an  orthonormal  basis  for 
Vj  and  is  an  orthonormal  basis  for  Wj, 

Then  L^(R)  =  that  is,  for  any  function 

/  E  L^(R),  if  we  let  Pj  denote  the  orthogonal  projection 
Z/^(R)  V},  then  Pj  converges  to  /  in  the  norm. 

lim  Pjf  =  /  (14) 

J— ►oo 

The  Mallat  identity  asserts  that  ^7 

Vj  =  Vo  0  Wo  0  *  ■  -  0  Wj-i  (15) 

If  one  assumes  an  expansion  of  the  form 

/i(x)  =  ^/>.Vi(x)  eVj  (16) 


then  one  can  also  write 

= E  + E  E  /‘’VK*) 

iez  .=0  itez 
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where 


In  '8  a  theorem  about  the  wavelet  interpolation  of 
a  function  in  2-D  using  only  the  scaling  function  is 
proven. 

Theorem  i4ssume  the  function  f  g  C^(fl),  where 
ft  is  a  bounded  open  set  in  R^.  Let  for  j  g  N, 

P.»€A  ^  ' 

where  x,y  g  ft,  the  index  set  A  =  {t  g  Z  :  supp{4>\)  n 
ft  0}  and  c  :=  /  x4>dx  =  ^  kh^.  Then 

\\f  -  r\\LH^)<C{llVf 


\\f  -  f%nm  <  c/v 

where  C  is  a  constant  depending  only  on  the  diameter 
of  rtj  N,  and  the  maximum  modulus  of  the  first  and 
second  order  derivatives  of  f  on  fl. 

In  any  application  involving  the  approximate  solu¬ 
tion  of  partial  differential  equations  using  wavelets,  it  is 
essential  to  derive  accurate  representations  of  the  weak 
form  of  the  governing  equations  in  the  wavelet  bases. 
For  some  classes  of  spline  wavelets,  the  derivation  of 
differentiation  formulae  is  rather  straightforward.  How¬ 
ever,  the  smoothness  properties  of  Daubechies  wavelets 
are  rather  abstract,  the  calculation  of  their  derivatives 
is  somewhat  delicate.  This  task  has  been  studied  in 
^9,^0,^  1  for  exact  quadratures.  The  effect  of  trunca¬ 
tion  error  in  the  differentiation  of  Daubechies  wavelets 
is  discussed  in  ^2 

The  approach  taken  in  this  paper  employs  transla¬ 
tions  and  dilations  of  the  scaling  function  only  as  in  ^8, 
as  opposed  to  using  the  associated  wavelet  bases  them¬ 
selves  ^9.  As  in  the  last  section,  if  a  function  F{x)  has 
the  approximate  representation  F^{z), 

f^(x)  =  EW(x)  (20) 


Then  one  can  differentiate  this  expression  to  achieve 

(21 


fk^  ^  0  <  t  <  7 

1  '  .  (18) 

E  0  <  t  <  j 

Conversely,  if  one  knows  the  coefficients  of  the  expan¬ 
sion  at  level  t,  the  coefficients  at  level  t  —  1  are 

^  ^(//  +  //“^’^ajb-2/),  0  <  t  <  7 


Consequently,  the  solution  for  the  derivative  of  F^{x) 
reduces  to  the  calculation  of 


for  all  p  E  iV  (22) 


However,  one  can  always  expand  the  derivative  in  terms 
of  the  basis  functions 


(<^^(3;))  =  Ylo:k,p<l>i{x) 


(23) 


Taking  the  L^-inncr  product  of  the  above  equation  with 
an  arbitrary  basis  function  (}>l{x) 

«r.p  =  J  <Prix)  dx  (24) 

and  by  using  a  change  of  variables  one  can  show  that 


One  can  derive  the  weak  formulation  of  the  plate 
equations  using  a  domain  embedding  technique  by  first 
defining  the  cost  functional  Jq  by 

Mv)  =  \ei  I  0(x)  p  dx  +  ifc,  I  J,(x)  r  dx 


Oir,p  =  V  J  +  r -p)d(,  ,  ^  =  2^x+p 

This  equation  Is  one  example  of  the  general  two-term 
coefficients  defined  by  ^0 

-  ^2))  dx 

In  particular, 

Q^r.p  =  2-^Aq’p_^  (26) 

so  that  the  final  expression  for  the  derivative  of  F^{x) 
becomes 


P6A  r  ^27) 

p€A  r 

For  the  calculation  of  numerical  boundary  measure 
explained  in  later  section,  we  also  need  the  three-term 
connection  coefficients  defined  by  [Latto] 

^h.b./3  “ 

~  ~  -  ^2)  d^ 

(28) 

It  is  important  to  note  that  these  connection  coeffi¬ 
cients  can  be  calculated  once,  stored  in  a  program  and 
used  repeatedly.  Detailed  discussions  of  their  calcula¬ 
tion  can  be  found  in  ^9,^  0,^  1. 


(IV)  Domain  Embedding,  Wavelet  Galerkin  Methods 

The  domain  embedding  technique  has  been  em¬ 
ployed  extensively  in  computational  mechanics  ,  and  in 
particular,  in  studying  subdomain  decomposition  meth¬ 
ods  .  In  this  section,  the  wavelet  galerkin  penalty  for¬ 
mulation  for  the  one  dimensional  beam  equation  is  de¬ 
rived.  The  extension  to  include  aerodynamic  terms  in 
the  panel  flutter  model  and  the  associated  equations 
for  the  control  design  of  flutter  suppresion  of  plates  are 
subsequently  straightforward  to  derive.  Consequently, 
the  extension  of  the  wavelet  galerkin  approximation  of 
the  beam  equation  incorporating  aerodynamic  terms  to 
model  panel  flutter  are  simply  stated  without  proof. 
The  interested  reader  is  referred  to  ^2  for  the  details. 


“  /  y{^)f{x)dx 
JD 


To  enforce  the  boundary  conditions  it  is  necessary  to 
define  the  trace  operator  7  that  evaluates  the  displace¬ 


ment 


and  the  slopes  at  the  boundary 


By  introducing  the  penalty  terms  on  the  boundary 


lo{y) 


a  modified  functional  can  be  defined 


(32) 


—  *^o{y)  d"  Jciy) 

The  Euler  equations  for  the  modified  functional  J{y) 
define  the  necessary  equations  for  equilibrium  of  the 
plate.  In  other  words,  by  taking  the  Gateaux  differen¬ 
tial  of  Jo{y) 


<  DJoiy),v  >=  El  0(x)0(x)dx 

+fci  j  y{x)v{x)dx~  I  v(x)f(x)dx  (34) 
Jd  Jd 

and  by  taking  the  Gateaux  differential  of  Je{y) 

<  DJ^{;y),v  >=  i  (7o.i(y)  -  y”)  7o,i(«) 

(7o,2(y)  -  y®)  7o.2(t<) 

(7i.i(y)  -yi)  71,1  (^) 

(71 ,2 (y)  -  yl)  71.2(1^)  (35) 


the  governing  equations  in  weak,  or  variational  form, 
become 


+  ki  J^y{x)v{x)di 


+-lo,i{y)lo,i{v)  +  ^7o,2(y)7o,2(«^) 
+-7i.i(y)7i,i(w)  +  -7i,2(y)7i,2(v) 


+  “ya7o,i(w)  +  7o.2(v) 

+jyi  7i,i(v)  +  ^yl  7i,2(v)  (36) 

The  galerkin  approximation  of  the  variationaJ  equations 
are  obtained  by  substituting  the  expression  for  the 
approximate  solution  and  trial  function 


v  =  '^vi<f>i{x) 


(37) 


With  this  substitution,  the  approximate  equations  be¬ 


come 


^  Jo 

+“Ej('  A-'W^T  yHt) 

E^  <^f(^)<P{(^)dxyf(t) 

+ci  ^  4>i  {x)4>\  (x)dx  yf  (t) 

•hki  ?/  d>i  ix)4>'l  {x)dx  yf  (t) 

■*■7  E  {  (a)  +  <f>f  ib)<l>f  (a)}  yf{t) 

i 

(»}»/(.) 

=  E  E  ^ir  [  d>i  (a;)^  {<pf  (a;))  dxui{t) 

i  p  Jo  ox 


+72^0  («)  +  1  y“  <j>f  (b) 


(38) 


To  obtain  the  final  equations  in  the  wavelet  basis,  one 
need  only  substitute  the  equations  derived  in  the  last 
expression  for  differentiating  wavelets  in  terms  of  the 
connection  coefficients,  or  wavelet  quadratures,  A^j 


E  (^b«-  y.- (0  +  (ci  -  6;,-  yf  (O) 

+  E  (h6„+2^-'EIAlf  -  A?;/)  yf(t) 

<f>f  {b)<i>f (a)}  yf{t) 
i 

+7  E  E  E  { <■“  A-:“  ^fc(«)<^r(«) }  yf{t) 

i  k  r 

+7EEE{  ^l:kA]::^<f>i{b)<f>fib)}  yf{t) 

i  k  r  ^ 


=  E  E  K^raihfr  «•■(<) 

t  r 

+“yS  («)  +  ^  y“  <l>f  (b) 

+ 7  E 1^'“  ^/.r  d>f{a)  +  ^yl  <6;!  (6)  I  (39) 

These  final  equations  have  several  properties  that 
make  them  attractive  from  a  computational  viewpoint: 

(i)  The  various  coefficient  matrices  are  highly  sparse, 

(ii)  the  wavelet  quadratures  are  calculated  once,  prior 
to  the  simulation,  and  stored  for  application  to  any 
partial  differential  equation, 

(iii)  rapidly  convergent  diagonal  preconditioning  meth¬ 
ods  can  be  derived  for  representations  of  differential 
operators  in  the  wavelet  basis  and 

(iv)  preliminary  study  by  the  authors  have  shown  that 
multigrid  wavelet  procedures  are  rapidly  convergent 


However,  it  should  be  noted  that  the  extension  of  the 
above  equations  to  more  general  domains  is  much  more 
difficult  than  appears  in  the  one  dimensional  case.  In 
the  case  of  arbitrary  domains  in  two  dimensions,  the 
penalty  term  that  must  be  evaluated  has  the  form 


-  J  uvdS  (40) 

where  fl  is  the  boundary  of  the  plate  and  7  is  the  trace 
operator  that  evaluates  a  function  defined  on  il  at  its 
boundary.  The  next  two  sections  summarize  methods 
for  calculating  these  penalty  terms  for  arbitrary  bound¬ 
aries. 


(V)  Geometric  Measure  Integrals 

The  first  method  is  introduced  by  Wells  ^8.  Because 
it  forms  the  basis  for  numerical  studies  and  the  new 
method  presented  in  the  next  section,  it  is  discussed  in 
some  detail  in  this  section.  Let  Q,  be  an  open  bounded 
set  in  where  boundary  is  defined  by 

dD,  =  {(^,1/)  €  R^  :  F{x^y)  =  constant}  (41) 

Then  the  unit  normal  vector  il  along  the  boundary  dO. 
can  be  written  as 

VF 

”  “  ivf  I 

By  extending  n  to  R^  “smoothly”,  we  have  the  follow¬ 
ing  formula  for  the  arclength  of  Oil 

L[dH)  =  /  ds  =  j  n  '  nds  ^  j  divii  dx  dy 

Jan  JdQ  JQ 

=  /  XndivTi  dxdy  =  -  /  •  n  dx  dy 

where,  xq  is  the  characteristic  function  of 

More  generally,  by  the  same  derivation,  one  has  for 
any  integrable  function  /  defined  on  517,  after  extending 
/  to  R^ 

f  fds  =  ^  [  /II51I11  (44) 

JdQ  7r2 

where  ]|5ni|  =  •  iidxdy. 

Now  one  can  develop  the  wavelet  approximation 

of  boundary  integral.  The  wavelet  expansion  of  the 

characteristic  function  xq  and  F{x,y)  at  a  level  j  is 

Xr2(a;,y)  = 

P.9 

=  (45) 

P,9 

Using  the  connection  coefficients 

^I’m  '=  [  4>\<Pmdx  (46) 

J  R 

it  is  possible  to  derive  the  expansion  coefficients  of  the 
derivatives.  The  partial  derivatives  in  terms  of  wavelet 
epansions  are 


p,9  ^  '  P^<1 

^  =  E(f)'  Wy) 

p,q  ^  ^  P.9 

dF^  V-  fdFV  . 

p,^  \  /  p,7 


For  large  j,  one  can  accurately  approximate  the  values 
of  ~  and  at  the  interior  points  by  their  correspond¬ 
ing  coefficients,  i.e. 


dF  .  .  ^:fdFy 


dF  .  fdFy 


|VF|(p,,)«2<  ((St')'  I  +(f|^')' 

\  \  \dx  J  \\dy  J 

^  \  ^  ^  P.9/  \  '^  ^  '  P.^ 

Letting  n=  (n^,n„), 


lx  (p,  9) «  -  I  (50) 


therefore  the  expansion  of  n^,  Uy  are. 


■i=E 


((f  ):j 

(«)i. 


‘f’pi^Wqiy) 


.  ^ I  vy  J 

~[  ■  A\>{x)4>q{y) 

And  the  numerical  boundary  measure  ||9fl||  is 

\\dQ.\\  =  Vxn-ndxdy  (5; 

||512||  =  (II^^IDp,,  <!>p{xWq{y)dx  dy 


Substitute  wavelet  expansion  of  Vxfjj  and  n,  take 
inner  product  with  <l>l{x)^{y),  and  use  the  three  term 
connection  coefficients 


=  /  4>p{x:)4>k{x)^s{x)dz 
J  R 


then  we  obtain 


[(li) 

(«);„] 

^P,*:,«^9.M 

(s)'.,.' 

r 

Now  wc  have  a  formula  for  calculating  the  arclength 

Lidn)  =  -  f  E  (ll^^^ll)L 


(56) 


(VI)  Positive,  Symmetric  Boundary  Mea¬ 
sures 

As  will  be  seen  in  Section  (VII),  while  the  numeri¬ 
cal  boudnary  measure  described  in  the  last  section  can 
be  effective  in  direct  calculation,  numerical  examples 
have  shown  that  it  docs  not  preserve  the  positive,  sym¬ 
metric  form  of  the  terms  arising  penalty  contribution 
to  the  variational  equations.  Specifically,  for  boundary 
integrals  having  the  form 


uvdS 


(57) 


it  is  possible  to  derive  a  numerical  scheme  that  retains 
the  positive,  symmetric  attributes  of  the  penalty  term. 
By  formula  (??),  one  can  always  write 


f  Vxfi  •  ndxdy 


(58) 


when  /  is  extended  to  all  of  R^.  To  obtain  a  symmetric, 
positive  measure,  one  approximates  the  boundary 


=  {{x,y)  e  =  l}  (59) 


and  the  normal  derivative  to  the  approximate  boundary 


n  ^  ft^  = 


11  11 


(60) 


By  assuming  the  form  of  the  solution  and  trial  functions 
to  be  a  linear  superposition  of  wavelet  scaling  functions. 


(61) 


the  final  form  of  discrete  approximation  for  domain 
embedding  is  obtmned 


^ih^id 


(62) 


the  four  term  connection  coefScients  defined 

to  be 

K'li'^d=  f  4>i{x)<i>iix)<j>i{x)4,-’j{x)dxdy  (63) 

The  calculation  of  the  four  term  connection  coefficients 
is  discussed  in  ^3. 


(VII)  Numerical  Examples 


Several  numerical  examples  have  been  carried  out 
to  verify  the  accuracy  and  consistency  of  the  formula¬ 
tion  for  the  1-D  wavelet  galerkin  approximation  of  the 
panel  represented  in  eqs.  (39).  Figures  (1)  through 

(6)  depict  the  solution  of  the  boundary  value  problem 
without  aerodynamic  terms  to  study  the  overall  conver¬ 
gence  rate  of  the  wavelet-galerkin,  embedded  domain 
problem.  In  the  first  simulation,  shown  in  figures  (1) 
and  (2),  a  clamped  panel  has  been  modeled.  Figure  (1) 
depicts  the  solution  on  the  full  domain,  while  Figure 
(2)  illustrates  the  solution  on  the  physical  domain  for 
several  resolution  levels.  Similar  results  are  shown  in 
Figures  (3)  and  (4)  for  a  simply  supported  panel,  while 
Figures  (5)  and  (6)  depict  the  results  for  a  clamped 
plate  subject  to  a  constant  transverse  load.  Figures 

(7) ,  (8)  and  (9)  show  the  asymptotic  rate  of  conver¬ 
gence  plotted  as  a  function  of  the  resolution  level.  As 
expected,  the  convergence  rate  is  0(2'^). 

Inasmuch  as  the  numerical  boundary  measure  is  de¬ 
fined  in  a  rather  abstract  manner.  Figures  (10)  through 
(12)  graphically  illustrate  the  form  of  the  boundary 
measure  and  its  constituents.  For  example.  Figure  (10) 
illustrates  the  characteristic  function  of  the  circle(a:^  -h 

=  100)  at  the  sampling  level  j=l.  Figure  (4)  de¬ 
picts  the  boundary  measure  |1  itself.  Clearly,  all  the 
constituent  functions  comprising  the  boundary  measure 
are  highly  sparse,  as  is  expected  in  wavelet  calculation 
of  localized  phenomena.  The  error  in  calculating  the 
arclength  of  the  perimeter  of  the  circle  is  depicted  in 
Table  (1)  and  agree  well  with  results  presented  in  ^8. 
Figure  (11)  is  a  characteristic  function  of  the  airfoil  and 
the  numerical  boundary  measure  ||9fll|  is  shown  in  Fig¬ 
ure  (12).  From  these  boundary  measures,  the  approx¬ 
imated  arclength  of  the  boundary  of  the  airfoil  for  the 
level  j  =  4  is 

L  =  8.219033 


Future  work  will  include  the  detailed  study  of  the 
coupled  fiuid-structure-control  interaction  problem  as¬ 
sociated  with  the  control  of  panel  flutter  using  a  wavelet 
galerkin  formulation.  In  addition,  fast  multigrid  algo¬ 
rithms  for  the  solution  of  the  resultant  equations  of  mo¬ 
tion  have  been  developed. 


Intheaboveequation^ 
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Figure  (2) :  Embedded  Domain, 
Clamped  Plate  Solution 


Figure  (4) :  Embedded  Domain, 
Simpy  Supported  Plate  Solution 


Figure  (6) :  Embedded  Domain, 
Clamped  Plate  Solution 
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Figure  (7)  Asymptotic  Convergence  Rate,  Clamped  Plate  1 
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Figure  (8)  Asymptotic  Convergence  Rate,  Hinged  Plate  1 
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Figure  (9)  Asymptotic  Convergence  Rate,  Clamped  Plate  2 


Figure  (12)  Characteristic  Function,  Airfoil  Figure  (13)  Boundary  Measure,  Airfoil 


Table  (1)  Approximated  Arclength  for  Circle  =  IQO 


Level 

Exact  Length 

Approximated  Length 

Error 

i  =  0 

62.8319 

62.812821 

-0.030295  % 

j  -1 

62.8319 

62.815614 

-0.0258435  % 

j  =  2 

62.8319 

62.809639 

-0.0353548  % 

i  =  3 

62.8319 

62.842482 

+0.0169165  % 

j  =  4 

62.8319 

62.831550 

-0.0005570  % 

Attachment(2) 
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ABSTRACT 

Tlie  investigators  present  a  study  of  dynamic  and 
aeroelastic  response  of  structures  which  evolve  due  to 
damage.  Aeroelastic  response  is  shown  to  be 
dependent  upon  the  distribution  and  accumulation  of 
damage.  In  turn,  the  damage  is  dependent  upon  the 
presence  of  the  aerodynamic  loads.  Dynamic 
characteristics  are  unique  to  the  coupled  damage  / 
aeroelastic  system  and  are  developed  as  part  of  the 
solution  methodology.  In  this  study,  the  damage  is  due 
to  the  natural  progression  of  microcracking  of  the 
composite  structure;  yet,  the  control  model  presented 
is  appropriate  for  distributed  actuation  systems.  Tlie 
stability  boundary  for  aeroelastic  flutter  and 
divergence  evolves  due  to  damage.  Control  design 
based  upon  the  min-max  control  theory  is  presented 
which  addresses  model  uncertainties. 


INTRODUCTION 

Aeroelasticity  is  defined  as  the  phenomenon 
resulting  from  the  interaction  between  aerodynamic, 
inertial,  and  structural  forces.  A  distinguishing 
characteristic  of  aeroelasticity  is  the  presence  of  time 
dependent  and  displacement  dependent  aerodynamic 
loads.  Traditionally,  the  elastic  forces  of  the  structure 
are  assumed  linear  and  the  stiffness  of  the  structure  is 
assumed  to  remain  constant  for  the  life  of  the 
structure.  However,  these  assumptions  are  not  valid  for 
structures  which  incorporate  composite  or  active 
materials.  Of  particular  interest  in  our  research  is  tlie 
change  in  the  aeroelastic  stability  due  to  the  evolving 
structure.  This  evolution  is  the  result  of  natural 
degradation  of  the  structure  due  to  a 
naturally-occurring  progression  of  damage.  Tliis 
damage  accumulation  is  dependent  upon  the  history  of 
the  aerodynamic  loads  and  the  resulting  defonnations; 
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yet,  these  aerodynamic  loads  are  dependent  upon  the 
increased  flexibility  of  the  structure. 

The  investigators  examine  the  effect  of  damage 
accumulation  in  aeroelastic  structures  and  elucidates 
the  unique  behavior  of  the  dynamics  for  the 
evolutionary  structure.  The  primary  goal  of  the 
research  includes  the  development  of  the 
fluid-structure  model  with  structural  evolution.  Tlie 
model  pennits  the  examination  of  dynamic  response 
and  aeroelastic  stability  characteristics  of  the  evolving 
system,  as  well  as  damage  detection  and  system 
identification. 

Tlie  concept  of  actively  suppressing  dynamic  or 
aeroelastic  response  by  using  smart  materials  is  being 
pursued  by  several  researchers.  Control  theory  is 
required  which  is  insensitive  to  partial  loss  of 
distributed  actuation  systems,  classes  of  unmodelled 
damage  effects,  evolution  of  the  material  constitutive 
laws,  or  unmodelled  nonlinear  aerodynamic 
influences.  Hence,  another  goal  is  the  development  of 
control  theory  for  the  evolutionary  structure;  these 
activities  are  addressed  in  a  related  paper. 

A  method  has  been  developed  to  treat  the 
interaction  of  aerodynamic  forces,  dynamics  of  the 
structure,  and  damage  accumulation.  Presently,  the 
method  is  developed  for  the  specific  case  of  panel 
flutter.  Panel  flutter  is  defined  as  self-excited 
oscillatory  motion  caused  when  panels  are  exposed  to 
high-speed  airflow  on  one  side  of  the  structure  (see 
Figure  1).**^  Few  researchers  have  examined  the 
effects  of  evolution  on  aeroelastic  phenomena. 

Dowell  briefly  described  the  change  in  the 
aeroelastic  stability  arising  from  the  degradation  of  the 
isotropic  panel  structure.^  Dowell  suggested  that 
nonlinear  flutter  analysis  would  permit  the  prediction 
of  the  fatigue  life  of  the  associated  panel  structure. 
Dowell  further  suggested  that  conventional  fatigue 
data  could  be  used  for  estimating  the  fatigue  life  of  the 
panel  structure  specifically  at  flutter  conditions.  Xue, 
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Tlie  governing  equations  for  the  linear  aeroelastic 
behavior  of  the  plate  become 

ph[M]{w}+g*[CA]{w}+[[K]+;^*[K;i]]{w}  =0 

where  p  is  the  density  of  the  material,  and  h  is  the 
plate  thickness.  We  also  introduce  the  parameters 
X*  =  —===-2Lnd  O'*  =  ^  for  convenience. 

Tlie  matrices  \_Ca\  and[-^^J  represent  unsteady 
aerodynamic  contributions. 

We  assume  a  solution  for  these  equations  of  the 
fonn  {w}  =  where  Q  =  5  +  /a).  The  equations 

are  written  in  nondimensional  form  as 

[[K]  +X[Ka]-k[M]](w}  =  0 

where 

^  =  =  plate  length,  Dr  -  reference  stiffness, 

Ka  =  K  =  K.  The  eigenvalue  is  defined  as 

K  = -(^)  ' .  where  a.  =  =  ^. 

and  ]i  =  for  the  case  M  »  1 .  Tlie  eigenvalues 

are  complex  and  dependent  upon  flow  conditions.  The 
matrix  \^K  includes  the  effect  of  damage  and  is 
dependent  upon  the  prior  load  history. 

Damage  Model  for  the  Structure 

Composite  materials  are  used  for  aircraft  and 
spacecraft  applications.  Tlie  ability  for  composite 
materials  to  provide  lightweight,  response-tailored 
structures  has  been  well  documented.  Adversely, 
microstructural  damage  is  both  inherent  and 
unavoidable  in  these  materials.  Sources  of  damage 
include  matrix  cracking,  delamination,  fiber  breakage, 
or  fiber-matrix  debonding.  As  a  consequence,  the 
research  to  characterize  this  damage  has  been  quite 
active  resulting  in  several  analytically  and 
experimentally  based  damage  models. 

We  use  a  continuum  mechanic  approach  for  the 
damage  model.  Tliis  method  follows  the  general 
formulation  given  by  Coleman  and  Gurtin  as  described 
in  their  work  on  thermodynamics  with  internal  state 
variables  (ISV).’  Tlie  method  proceeds  from  the 
assumption  that  IS  Vs  can  be  implemented  to  describe 
the  state  of  damage.  Damage  effects  are  introduced  at 
the  ply  level;  thus,  the  ISV  method  is  independent  of 
ply  stacking  sequence.  An  important  characteristic  of 


the  ISV  approach  is  that  an  elastic  medium  with 
damage  is  treated  as  an  equivalent  homogeneous 
material. 

The  damage-dependent  stiffness  components 
{  described  in  terms  of 

the  damage-dependent,  angle-ply  stiffness,  [q"].* 
Talerja  defines  a  parameter,  to  describe  damage 
within  [q^].  The  stiffness  of  the  composite  laminate 
with  damage  is  found  as 

,  B' J  ,  D' j)  =  j  (1 ,  Z,  z2)  dz 

where 

Q'a  =  Q'lj  ,  E2(.Q ,  Vi2(q) ,  Gi2©] 

The  properties  for  the  principal  directions  of  a 
lamina  are  given  by 

Bl  =  £?  +  2^(C3  +  C8(v?2)2  -  Ci6  v?2) 

E2=£?  +  2^(C8  +  C3(v?J  -Ciev^i: 

V12=V?2  +^(-iZ^)(Ci6-2C8V?2) 

<^12  =^12  +  2^Ci3 

where  £|,  ^12^  ^12  properties  of  the 

undamaged  lamina.  C3,C8,Ci3,  and  C]6  are  constants 
phenomenological  in  nature.  Tlie  damage  parameter^ 
represents  the  crack  density,  the  magnitude  of  which  is 
subject  to  a  damage  growth  law. 

The  damage  growth  model  incorporates  the  efforts 
of  Allen,  et  al  and  Talreja.^'*^ 

6Cm/i  =  T,  CX^) 

where  is  the  set  of  p  IS  Vs,  ejy  is  the  strain,  T  is 
the  temperature,  and  the  overdot  represents  the 
derivative  with  respect  to  time. 

The  ISVs  are  defined  as 

OLy  =  ^  j  UifljdS 
^  s 

where  a,y  are  components  of  the  ISVs,  Vl  is  the  local 
volume  in  which  homogeneity  is  assumed;  Ui  is  the 
crack-face  displacement,  and  rij  is  the  normal. 


Fig.3  The  evolution  of  stiffiiess  for  the  aeroelastic 
structure. 

use  of  adaptive  control  strategies.  Adaptive  control  is 
imperative  for  models  such  as  those  described  herein. 
The  basis  for  this  control  strategy  will  be  an  extension 
of  the  adaptive  Hamiltonian  control  theory. 

Although  some  efforts  have  been  made  to  derive 
active  control  strategies  for  aeroelastic  structures  using 
smart  materials,  this  research  is  in  the  early  stage  of 
development.  The  authors  are  aware  of 


Fig.  4.  The  natural  frequencies  of  a  composite  plate 
with  the  progression  of  microstructural  damage. 


no  research  that  treats  the  fundamental  issue  of 
robustness  of  control  desigii  witli  respect  to  actuator 
degradation  expected  during  routine  operation,  let 
alone  the  possibility  of  sensor  or  actuator  failure. 
While  numerous  research  papers  on  tlie  use  of  smart 
materials  for  active  control  are  available,  the  works  of 
Scott  and  Rogers  arc  particularly  pertinent. In 
Scott,  an  LQR  strategy  is  employed  to  design 
piezoelectric  controls  to  reduce  pane!  flutter  and  raise 
the  value  of  dynamic  pressure  associated  with  flutter. 
The  analytical  results  are  extremely  promising.  On  tlie 
other  hand,  the  work  of  Rogers  uses  shape  memory 
alloys  to  reduce  crack  growth  in  composites.  One 
potential  drawback  with  the  approach  taken  by  Scott  is 
that  the  LQR  design  strategy  docs  not  account  apriori 
for  injected  disturbances  or  material  change  which  will 
develop  in  fatigue  or  damage.  Likewise,  Rogers 
clearly  demonstrates  that  shape  memory  alloys  have 
the  potential  to  significantly  reduce  local  stress 
concentrations;  however,  an  active  control  strategy 
based  upon  some  optimality  principle  does  not  exist. 
Furthenuore,  the  approach  in  Rogers'  work  deals  with 
degradation  of  the  composite  lajuinate  while  the 
integrity  of  the  control  is  unaffected. 

The  design  of  control  strategies  for  the  aeroelastic 
model  with  damage  accumulation  cau  be  interpreted  as 
a  two  player  game.  The  simplest  explanation  is 
achieved  by  examining  the  prior  governing  equations 
but  including  control.  In  this  case,  the  two  right  hand 
side  tenns  have  "adversarial"  roles, 

[M]{w}  +[K]{w}  =[[Kd]  +  [A]]{w}  +[B]{u} 

The  stiffness  damage  and  aeroelastic  tenns, 
C[^d]+[A]]  ,  tend  to  destabilize  the  system,  while  the 
control  terms,  [B],  are  selected  to  offset  this  effect. 
The  destabilizing  terms  represent  localized  damage 
from  microcracking  or  delamination,  actuation  loss, 
unmodelled  disturbances,  or  umnodelled  nonlinear 
effects. 

The  investigators  have  extended  the  above 
rudimentary  model  to  account  for  distributed 
parameter,  or  infinite  dimensional,  models  of  system 
dynamics  for  classes  of  disturbances  and  uncertainty. 
A  more  detailed  discussion  of  this  approach  can  be 
found  in  related  papers.'®  ’’*^® 
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Fig,  6.  The  mode  shapes  of  the  cantilevered  composite  plate.  Some  modes  are  influenced  by  aerodynamic  loads. 
The  primary  bending  and  twist  mode  coalesce  to  fonn  the  aeroelastic  modeionfirst 


Of  particular  interest  is  the  reduction  in  the  critical 
dynamic  pressure  for  the  case  of  damage  growth  due 
to  the  stress  field,  but  the  increase  in  the  critical 
dynamic  pressure  for  the  case  of  uniform  damage. 

The  mode  shapes  of  the  aeroelastic  structure  are 
dependent  upon  the  dynamic  pressure.  The  first  tliree 
mode  shapes  are  shown  in  figure  6  for  increasing  flow 
conditions.  At  flutter,  the  first  two  modes  coalesce  to 
form  a  single  aeroelastic  mode.  Tlie  third  mode  is  not 
influenced  at  these  conditions.  Although  the 
frequencies  show  a  significant  change  due  to  damage, 
the  differences  in  the  mode  shapes  due  to  damage  are 
not  remarkable. 

The  aeroelastic  boundary  is  dependent  upon 
damage  and  fiber  orientation.  The  aeroelastic  stability 
boundary  is  presented  in  fig.  7  for  several  orientations. 
Unifonn  damage  is  assumed  to  illustrate  the  effect  of 
damage  on  this  boundary.  Tlie  minimum  boundary  is 
presented:  the  structure  is  limited  by  flutter  for  certain 
fiber  orientations  fiber  angles;  the  structure  is  limited 
by  divergence  for  other  fiber  orientations.  Damage  has 
a  significant  effect  on  the  stability  boundary. 


Damage  is  not  unifonnly  distributed;  rather,  the 
distribution  is  dependent  upon  the  stress  field  which  is 
not  symmetric  due  to  the  presence  of  aerodynamic 
loads.  Tlie  stress  distribution  for  the  composite  lamina 


Fig.  7.  The  aeroelastic  boundary  for  the  composite 
structure.  Flutter  and  divergence  depend  upon  the 
lamina  orientation  and  microstructural  damage. 
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Abstract 


This  paper  presents  a  novel  strategy  for  obtaining 
finite  dimensional  approximations  of  Riccati  operators 
arising  in  game  theoretic  control  problems  on  Hilbert 
spaces.  While  extensive  progress  has  been  made  in  the 
utilization  of  the  newly  developed  field  of  wavelet  anal¬ 
ysis  in  signal  and  image  processing,  for  less  progress  has 
been  made  in  applicationa  to  computational  mechanics. 
This  fact  is  surprising  considering  the  inherent  multires¬ 
olution  properties  of  wavelet  analysis.  This  paper  de¬ 
rives  a  wavelet  galerkin  formulation  of  a  specific  class 
of  computational  control  problems:  zero  sum  dynamic 
games  on  Hilbert  spaces.  This  work  extends  previous 
research  in  that  not  only  are  convergence,  exponential 
stability  and  robust  stability  of  soft  constrained  differ- 
tial  games  on  Hilbert  spaces  considered,  but  quantita¬ 
tive  numerical  conditioning  results  are  also  presented. 
In  fact,  a  precise  bound  is  given  on  the  error  induced  by 
truncating  wavelet  quadratures  used  in  the  weak  formu¬ 
lation  of  computational  control  problems.  The  theory 
presented  is  verified  by  modeling  a  robust  control  strat¬ 
egy  for  a  beam  subjected  to  piezoceramic  actuation  and 
structured  uncertainty. 


N  omenclat  ure 


{ho, h2W-i}  Daubechies  coefficients  set  for 
scaling  function 
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system 
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€i  Structured  perturbation  scaling 

Si  Structured  perturbation 

Introduction 

With  the  ground  breaking  work  of wavelets  and 
multiresolution  analysis  have  become  a  predominant 
area  of  research  in  scientific  computing.  To  a  large  ex¬ 
tent,  most  work  to  date  has  focussed  on  image  process¬ 
ing  and  signal  processing^.  StiU  the  potential  for 
major  advances  in  computational  techniques  for  par¬ 
tial  differential  equations^  and  integral  equations^  have 
also  been  noted.  Just  ets  recently,  the  advantages  of 
multiresolution  analysis  in  certain  applications  in  com¬ 
putational  mechanics  has  been  established®’^. 

In  the  computational  mechanics  area,  the  finite  ele¬ 
ment  method  has  become  the  standard  by  which  other 
numerical  techniques  are  measured.  One  general  fea¬ 
ture  of  the  finite  element  method  is  that  as  the  mesh  size 
h  becomes  smaller,  the  resulting  matrix  of  the  algebraic 
equations  has  a  condition  number  increasing  at  the  rate 
of  (^) .  This  effect  can  result  in  very  slow  convergence 
rates  when  classes  of  iterative  techniques  are  used  to 
solve  the  resulting  algebraic  equations.  On  the  other 
hand,  and  have  shown  that  the  wavelet-based 
partial  differential  equation  formulations  yield  a  ma¬ 
trix  which  is  proven  to  have  a  nearly  constant  condition 
number  with  simple  diagonal  preconditioning,  indepen¬ 
dent  from  the  mesh  size.  When  one  considers  its  rapid 
convergence  properties  and  the  relatively  simple  com¬ 
putational  implementation®,  it  is  apparent  that  there 
is  an  enormous  potential  for  significant  development  of 
wavelet-Galerkin  methods  in  computational  mechanics. 

In  order  to  solve  differential  equations  of  any  type 
using  a  Galerldn  method,  one  eventually  needs  the  ap¬ 
proximation  of  the  derivative  of  some  basis  functions.  If 
one  uses  a  wavelet  basis,  such  as  the  Daubechies  wavelet 
system,  this  requires  the  calculation  of  connection  co¬ 
efficients  having  the  form 

Ajb  =  ^(x)<^(a:  -  k)  dx 

where  ^  is  a  wavelet  scaling  function.  While  the  connec¬ 
tion  coefficients  for  some  wavelet  systems  can  be  calcu¬ 
lated  in  closed  form  or  to  any  prescribed  accuracy 

using  numericeil  techniques  truncated  quadratures 
must  be  used  in  practice. 


Wavelet  Interpolation 

In  this  paper  Daubechies  wavelets  are  employed 
exclusively,  which  are  introduced  in  Let  <j)  be  the 
Daubechies  scaling  function  of  order  in  ^),  and 


let  {ho,  hi,  •  *  • ,  h2N-i}  be  the  corresponding  coefficients 
for  the  2-scale  difference  equation 

2^'-! 

(f>{x)  =  y/2  ^  hfc<;6(2i  -  h),  <^(x)  G  L^(R)(2.l) 

fc=0 

The  coefficients  {ho,  hi,  •  •  • ,  h2Ar-i}  satisfy 
2N^l 

V2  hk  =  2  .  (2.2) 

kzzO 

2N-1 

=  ^0,/,  I  E  Z 

jb=0 

where  Si^rn  is  the  Kronecker  delta  function.  The  coeffi¬ 
cients  {hk}  for  various  N  are  given  in  By  choosing 

bk  :=  (~l)^'^^h2Ar_i_fc  (2.3) 

[Daubechies]  derives  a  compactly  supported  wavelet 
function  from  the  scaling  function  ^(x)  via 

2N-1 

V’(*)  =  bk4>{2x  -  k),  i}>{x)  e  A^(R-)  (2.4) 

kzzO 

Since  the  coefficients  h^’s  are  firdte,  the  functions  ^  and 
Tp  both  have  finite  support  in  the  interval  [0,2iV'~l],  and 
they  are  orthogonal  to  themselves  and  to  each  other. 

One  can  define  the  wavelet  system  associated  with 
<j>  by 

<p{{x)  ~  Vl^4>{2^x  -  k)  (2.5) 

■^(x)  :=  2^^^'0(2^x  —  h) 

where  h  G  Z,  j  G  Z+(the  nonnegative  integers).  Then 
the  above  functions  have  following  properties  : 

supp  (^(*),  = [^,  ^  y — - 

IKIUr)=1  (2-6) 

X)(^(®)  =  2^,  x€R 

I  x"'*^(x)  dx  =  0  ,  m  =  0,v-,iV— 1 
Jr 

Moreover  4>y^  £  C®(^)  where,  for  example,  a(2)  « 
0.55,  o:(3)  «  1,088,  a(4)  «  xl,618.  Due  to  the  orthogo¬ 
nality  of  ^(x)  and  dilated  and  scaled  wavelet 

system  ^(x)  and  ^-^so  are  orthonormal, 

f  4>k{^)4>m{^)^^-  ^3,1' (2.7) 
Jr 

'  Jr 


If  one  defines 


or,  conversely 


Vj  ;=clos(span(<^j^(x),  k  j  fixed))  (2.8) 

Wj  ~clos(span('0{(a:),  k  j  fixed)) 

one  obtains  a  multiresolution  analysis  having  the  fol¬ 
lowing  properties  : 

(i)  •  C  V„_i  C  V„  C  V„+i  C  •  •  • 

(ii)  closi,(U„ezV„)  =  L"(R) 

(iii)  n„gzV„  =  0  (2.9) 

(iv)  =  V„_i©W„_i 

(v)  ©„6zW,-  :=-.-©W_i©Wo©Wi©---  =  i;2(R) 

Let  Pj  denote  the  orthogonal  projection  i/^(R)  — *  Vj. 
Then  Pjf  converges  to  /  in  the  norm, 


lim  Pjf  -  f 


J— oo 


(2.10) 


Let  P  {x)  be  the  projection  of  f  onto  Vy,  i.e. 

fH=^)  =  Pjf,  e-Vj  (2.11) 

From  the  MRA(Multiresolution  analysis)  property  (iv), 
p  (x)  can  be  decomposed  into  two  functions  such  that 


/•’(*)  =  /’  +  ^(®) 


(2.12) 


where 


/"■-^(x)  6  V,-_1  and  g^-^(x)  G  Wy_i  (2.13) 

If  one  continues  the  decomposition  M  times 

/i(x)=/>-"(x)  +  ff’-^(*) 

+  gi-"'+i(*)  +  ...  +  pi-i(x)  (2.14) 

Since  one  can  express  p{x)  in  terms  of  the  scaling 
function,  the  above  decomposition  can  be  rewritten  as 

r{x)=Y,44i^) 

k 

= E  4" + E  (2.15) 

k  k 

=  E4"^4""(*) 

k  k 

+"*+E‘^rM"^(®) 

k 

As  a  result,  one  can  derive  the  relationships  between 
the  coefficients  set  c][,  and 


4  ^  =  ;^E4^'-2»: 
4~‘  =  ;^E4«>*-2* 


(2.16) 


4  —  ;(^E  ^4  +  E4  (217) 

Approximations  to  the  coefficients  c]^  to  express  the 
function  P{x)  can  be  calculated  using  the  following 
quadrature  formula® 

< = (^) 


where  c  is  some  constant  defined  in  the  following  theo¬ 
rem. 


Theorem  Assume  the  function  f  G  C^(n),  where  Q  is 
a  bounded  open  set  in  IR?,  Let  for  j  G  N, 


==|-  E  /  (^. 

(2.19) 


X,y  efl 


where  the  index  set  A  =  {t  G  Z  :  supp(<^)  C\Ct  0} 
and  c  :=  f  x<f)dx  =  ^  Then 


and 


lh-/’lL.«»sc(i/2')= 


where  C  is  a  constant  depending  only  on  the  diameter 
of  Q,  N,  and  the  maximum  modulus  of  the  first  and 
second  order  derivatives  of  f  on  Q. 


Wavelet  Differentiation 


A  critical  feature  and  advantage  of  employing  a 
wavelet-Galerkin  basis  in  applications  to  computational 
mechanics  is  that  explicit  formulae  can  be  derived  for 
classes  of  differential  operators  of  the  connection  coeffi¬ 
cients  discussed  earlier.  Techniques  to  accomplish  this 
representation  are  discussed  in  The  approach 

taken  in  this  paper  employs  translations  and  dilations 
of  the  scaling  function  only  as  in  ®,  as  opposed  to  using 
the  associated  wavelet  bases  themselves.  As  in  the  last 
section,  if  a  function  f(x)  has  the  approximate  repre¬ 
sentation  P(x)f 


=  E4'4(*) 


(3.1) 


then  one  can  differentiate  this  expression  to  achieve 

=  («<*)) 

p€A 

Consequently,  the  solution  for  the  derivative  of  P  (x) 
reduces  to  the  calculation  of 

foraU  pGN  (3.3) 

However,  one  can  always  expand  the  derivative  in  form 
of  the  basis  functions 

^  (^(=))  =  (3-4) 

k 

Taking  the  X^-inner  product  of  the  above  equation  with 
an  arbitrary  basis  function  ^(x),  one  finds  that 

=  (3.5) 

By  a  change  of  variables  one  can  show  that 

«r.p  =2^'  J  +  r-p)di,(  =  Vx  +  p(3.6) 

This  equation  is  one  example  of  the  general  two-term 
coefficients  defined  by 

Kt  =  /  iH-  -  h))  dx  (3.7) 

In  particular, 

<^r,p  =  2-^  Aq*p_^  (3-8) 

so  that  the  final  expression  for  the  derivative  of  p  (x) 

becomes 

=  (3.9) 

p€A  r 

pGA  r 

It  is  important  to  note  that  these  connection  coefHcients 
can  be  calculated  once,  stored  in  a  program  and  used 
repeatedly.  Detailed  discussions  of  their  calculation 
can  be  found  in  and  Symbolic  manipulation 

programs  to  calculate  these  coefficients  to  arbitrary 
order  have  been  developed  by  the  author. 

Game  Theoretic  Control 

As  noted  earlier,  an  enormous  literature  has  accu¬ 
mulated  over  the  past  ten  years  regarding  the  conver¬ 
gence  of  finite  dimensional  approximations  to  the  in¬ 
finite  dimensional  linear  quadratic  regulator  problem. 


This  problem  is  formulated  by  considering  the  abstract 
Cauchy  problem  on  a  Hilbert  space  H 

i{i)  =  Ax{i) -i- Bu{t)  (4.1) 

x(0)  =  xo- 

where  U  is  o,  Hilbert  space,  B  G  X(f7,  H)  and  A  is 
the  infinitesimal  generator  of  a  strongly  continuous 
semigroup  T(i)  on  H.  By  defining  the  cost  functional 

•^o(«)  =  ^  +  (4.2) 

where  C  €  L[H,  H)  and  is  a  positive  definite  selfad- 
joint  bounded  operator  on  U,  the  LQR  problem  is  to 
find  u°eU  =  L^{0,  oo;  U)  such  that 

Jo{u°)  =  inf  Jo{u)  (4.3) 

subject  the  equations  of  dynamics.  However,  the  for¬ 
mulation  in  equations  (4.1)  and  (4.2)  does' not  account 
for  any  exogeneous  disturbance,  or  uncertainty,  in  the 
model.  To  obtain  a  theory  that  makes  provision  for 
disturbance  attenuation  and  uncertainty  for  infinite  di¬ 
mensional  models,  the  authors  have  extended  cetain  re¬ 
sults  from  flnte  dimensional  /  game  theoretic  con¬ 
trol  and  approximation  theory  for  the  LQR  problem  in 
equations  (4.1)  and  (4.2) 

The  theory  presented  in  considers  the  ab¬ 

stract  Cauchy  problem 

x(t)  =  (A  -h  AA)x(t)  -}-  (B  -h  AB)u(<)  (4.4) 

x(0)  =  xo, 

where  AA  G  X(iT,  H)  and  AB  G  X(?7,  H)  are  bounded 
linear  operators  representing  the  structured  uncertainty 
in  the  model  and  the  control  system.  Kurdila  and 
Fabiano^®  extend  the  finite  dimensional  results  in  by 
assuming  that  one  can  factor  the  perturbations  AA  and 
AB  as 

AA  =  DLaE  (4.5) 

AB  =  FLbG 

where  D,  X?  G  L{E,  H),  F  G  L{U,  H),  G  G  L{U,  U)  are 
known,  and  La  G  L{H,  H)^  Lb  E  L{U,  U)  are  unknown. 
The  robust  stability  result  in  and  can  be  derived 
by  introducing  the  fictitious  disturbances 

=  La  yi  (4.6) 

=  Lb  y2 

and  associated  observations 

2/1  =  Ex 

1/2  ==  Gu.  (4.7) 


These  equations  can  be  expressed  more  compactly  by 
letting  w  —  (it;i  ^v)2)  ^  H  X  H  and  y  =  (t/i  ,  ^2)  G  JT  X  17. 
In  this  case 


La  0 


0  Lb 


y 


(4.8) 


From  these  definitions,  the  system  model  with  struc¬ 
tured  uncertainty  can  be  interpreted  as  a  form  of  the 
original  system  with  injected  disturbance  and  can  be 
written  as 


x(t)  =  Ax{i)  -f  Bu(i)  -h  rto(t)  (^•^) 

a;(0)  =  xq. 


Here  F  G  L{H  x  HfH)  is  defined  by  = 

Dwi  +  Fit;2.  As  noted  in  (4.9)  can  be  treated  in 
the  differential  game  framework  by  considering  the  cost 

Ji{u,w)  =  ||Cx|^  + 

+|yp  -4h'  (4-10) 

KXtr  ^HxH  \  ^  ' 

/•OO 

— I 

y2  *  ‘vr  J 

where  Q  =  C*CaE*E,  R  ==  R+G*G,  and  VF  =  TTx^. 
With  this  soft-constrained  linear  quadratic  dynamic 
game,  one  can  associate  the  Hilbert  space  counterpart 
of  the  well  known  finite  dimensional  Riccati  equation 

A*U  4-  HA  -h  Q  -  n{BR^^B*  -  =  0.  (4.11) 

Reference  derives  the  following  result  extending  to 
infinite  dimensional  systems  the  robust  stability  theo¬ 
rem  proven  in 

Theorem  2.2  Suppose  that  (C,  A)  is  detectable.  If  H 
is  a  nonnegative  definite  solution  of  (4^1),  then  the 
feedback  control  u(t)  =  —R~^B*Jlx\t)  stabilizes  (4.4) 
for  all  IL^I  <7/  \Lb\  <  K- 

As  pointed  out  in  the  result  in^^  requires  that  (Q,  A) 
be  detectable,  while  this  result  only  requires  that  (C,  A) 
be  detectable. 

In  some  applications,  such  as  the  one  considered 
in  and  in  the  next  section,  standard  results  on  the 
convergence  of  approximations  to  the  infinite  dimen¬ 
sional  LQR  problem  [e.g.  can  be  applied  to  ap¬ 

proximate  the  game  theoretic  control  problem  discussed 
above.  Let  Le  a  family  of  finite  dimensional 

approximating  subspaces  of  H  such  that 


H=\jHn 

n=l 

and  suppose  that  An  G  L{En,Hn)t^  G  L{H,Hn),  and 
Qn  G  L{Hnj  TTn),  while  Pn  is  the  orthogonal  projection 


onto  Hn-  The  nth  approximating  algebraic  Riccati 
equation  can  be  written  as  an  operator  equation  on  Hn 

A;nn  +  AnTln  -  Un^n^r^Un  +  Qn  ==  0 
The  following  theorem  is  derived  in 

Theorem  2.3  Suppose  that  for  each  x  £  H  one  has 

(i)  Tn{t)PnX  ->  T{t)x  and  T;(t)PnX  T*(t)x 

(ii)  flnX  — ♦  Qx  and  Q*  Pn®  — ^ 

(in)  QnPn®  — ^  Q®  and  the  family  of  pairs  (An,  On) 
and  (AnjHn)  are  uniformly  stabilizable  and  uni¬ 
formly  detectable,  respectively  .Then  the  solution 
Hn  of  Lhe  finite  dimensional  approximation  al¬ 
gebraic  Riccati  equation  converges  strongly  to  H, 
the  unique  solution  of  the  game  theoretic  alge¬ 
braic  Riccati  equation 

lim  ||nx  —  HnPn®!!  — ^  0 

n — »oo 


Numerical  Examples 


As  a  numerical  example,  this  section  applies  the  the¬ 
ory  presented  in  the  paper  to  the  control  of  a  beam  with 
piezoceramic  actuators  and  subject  to  model  uncer¬ 
tainty.  In  practice,  the  model  uncertainty  may  be  due 
to  imprecise  actuator  characterization,  actuator  loss  of 
authority  due  to  material  degradation  or  damage.  Fig¬ 
ure  (3)  and  (4)  depict  a  redundantly  controlled  beam 
with  piezoceramic  actuators.  In  the  absence  of  model 
uncertainty,  the  governing  equations  characterizing  the 
dynamic  response  due  to  piezoelectric  actuation  can  be 
written 


N 


(5.1) 


where 

Ei{x)  =  JT(x  -  ai)  -  E{x  -  ft  )  (5.2) 

and  where  p  is  the  material  density,  h  is  the  laminae 
thickness,  El  is  the  bending  modulus,  ci  is  the  viscous 
damping  coefficient,  fci  is  the  stiffness  constant  of  the 
beam  foundation,  and  gi  is  the  piezoelectric  gain.  Each 
Qi  depends  upon  the  specific  laminate  geometry.  The 
interested  reader  should  see  or  for  details.  To 
account  for  several  regions  of  degradation  or  damage, 
structured  uncertcdnty  is  represented  in  the  system  as 


d^y  y^T^^y  ^y 

a?- + 


or,  alternatively,  as  injected  disturbances 


t  r 


4-  4-r  J-t 

t 

M 

+  ^ei7i(i)u;j(t),  (5.4) 

t 

The  utility  of  wavelet  approximation  for  control  prob¬ 
lems  becomes  evident  when  one  represents  the  control 
influence  functions  Hi{x)  and  uncertainty  distributions 
7i(x)  in  terms  of  their  wavelet  interpolates  as  defined 
in  section  (3) 

Hf{x)  =  J2Kr<l>r{-)  (5-5) 

r 

r 

By  substituting  the  galerkin  approximation 


t  r 

Results 

This  paper  has  presented  an  innovative  formulation 
using  a  multiresolution  wavelet  analysis  of  distributed 
control  problems  associated  with  soft  constrained  dif¬ 
ferential  games  on  Hilbert  space.  As  opposed  to  ap¬ 
proximate  optimal  controls  derived  using  linear  control 
theory,  the  game  theoretic  approximations  inherently 
provide  for  structured  uncertainty.  Numerical  verifica¬ 
tion  of  the  theory  has  demonstrated  that  the  formula¬ 
tion  is  well-suited  for  modeling  piezoceramic  actuation, 
and  suggests  avenues  of  research  for  the  control  synthe¬ 
sis  accounting  for  the  loss  of  actuation.  Future  research 
will  focus  on  the  conditioning  and  multigrid  solution 
of  the  multiresolution  analysis,  as  well  as  formulations 
of  fluid-structure  interaction  problems  with  robust  con¬ 
trol. 


/(==)  =  (5.6) 

r  References 


into  equation  (5.1),  multiplying  by  the  test  function 
V  =  {x)  ,  the  resulting  governing  equations  become 


<pa^)<l>i{^)dxyiii)  (5.7) 

E  ^  (®)  ^  W  (®))  yi  (*) 

+Ci  E  {x)dx  y(  (f) 

+*'1 E  (*)‘^®  yi 

=  E  E  Kr  <f>l  (®)^  (*))  dxUi{t) 


By  using  the  definition  of  the  connection  coefficients 
and  the  wavelet  differentiation  techniques  presented  in 
section  (3),  these  equations  reduce  to 


E  yi  (0 + Vi  (0)  (5-8) 
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ABSTRACT 

The  authors  present  a  study  of  the  dynamic  re¬ 
sponse  and  control  of  coupled  fluid-structure  systems 
in  which  the  primary  structure  or  distributed  control 
system  evolves.  These  systems  include  control  sensing 
and  control  actuation  elements  embedded  in  the  struc¬ 
ture.  Structural  systems  considered  arc  those  in  which 
damage  occurs  on  a  slow  time  scale,  or  damage  that 
occurs  due  to  microcracking  and  delamination.  Aeroe- 
lastic  response  is  shown  to  be  dependent  upon  the  dis¬ 
tribution  and  accumulation  of  damage,  which  itself  is 
dependent  upon  the  presence  of  the  aerodynamic  loads. 
Dynamic  characteristics  arc  unique  to  the  coupled  dam- 
age  /  aeroelastic  system  and  are  developed  as  part  of 
the  solution  methodology.  The  authors  establish  the 
importance  of  the  intrinsic  interdependence  of  stability 
and  damage,  and  establish  the  importance  of  deriving  a 
control  theory  for  uncertainties  in  the  evolving  system. 
Findings  indicate  significant  changes  occur  in  the  dy¬ 
namic  response  characteristics  due  to  damage.  These 
results  further  suggest  some  control  strategies  can  en¬ 
counter  difficulties  with  these  evolving  systems. 

INTRODUCTION 

Currently,  the  research  community  is  vigorously  ex¬ 
ploring  the  concept  of  actively  suppressing  dynamic  in¬ 
stabilities  by  using  ^^smart”  materials.  Furthermore, 
the  engineering  community  has  developed  new  materi- 

“  such  as  composite  materials,  piezoelectric  materi¬ 
als,  and  shape  memory  alloys  -  to  tailor  the  structural 
response-  The  primary  objectives  of  the  research  de¬ 
scribed  in  the  paper  include  the  development  and  verifi¬ 
cation  of  a  control  model  for  the  coupled  fluid-structure- 
control  system  with  struct  Ural /control  evolution.  Much 
of  the  related  research  does  not  consider  the  existence 
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of  damage  or  change  in  the  constitutive  behavior.  Con¬ 
trol  methodologies  arc  required  which  address  the  par¬ 
tial  loss  of  distributed  actuation  systems,  distributed 
damage  of  the  structure,  or  evolution  of  the  material 
constitutive  laws.  Control  theories  are  required  which 
are  fault  tolerant  with  respect  to  loss  of  either  sensing 
or  actuation  capabilities  over  portions  of  the  structure. 
The  investigators  examine  the  effect  of  damage  accumu¬ 
lation  in  aeroelastic  structures  and  elucidate  the  unique 
behavior  of  the  dynamics  for  the  evolutionary  structure. 
The  primary  goal  of  the  research  includes  the  develop¬ 
ment  of  the  fluid-structure  control  model  with  struc¬ 
tural  evolution.  The  model  permits  the  examination  of 
dynamic  response  and  aeroelastic  stability  characteris¬ 
tics  of  the  evolving  system,  as  well  as  damage  detection 
and  system  identification. 

Aeroelasticity  is  the  phenomenon  resulting  from  the 
interaction  between  aerodynamic,  inertial,  and  struc¬ 
tural  forces.  A  distinguishing  characteristic  of  aeroelas¬ 
ticity  is  the  presence  of  time  dependent  and  displace¬ 
ment  dependent  aerodynamic  loads.  Traditionally,  the 
elastic  forces  of  the  structure  are  assumed  linear  and 
the  stiffness  of  the  structure  is  assumed  to  remain  con¬ 
stant  for  the  life  of  the  structure.  However,  these  as- 
sumptions  are  not  valid  for  structures  which  incorporate 
composite  or  active  materials.  A  method  is  developed 
to  treat  the  interaction  of  aerodynamic  forces,  dynam¬ 
ics  of  the  structure,  damage  accumulation,  and  control. 
The  method  is  developed  for  the  case  of  panel  flut¬ 
ter  which  is  a  self-excited  oscillatory  instability  caused 
when  panels  are  exposed  to  high-speed  airflow  on  one 
side  of  the  structure^.  Few  researchers  have  examined 
the  effects  of  evolution  on  aeroelastic  phenomena,  Kim 
et.al^  and  Strganac  et.al^  present  the  development  of 
an  aeroelastic  model  in  which  damage  is  inherent. 

Dowell  briefly  described  the  change  in  the  aeroelas¬ 
tic  stability  arising  from  the  degradation  of  the  isotropic 
panel  structure^.  Dowell  suggested  that  nonlinear  flut¬ 
ter  analysis  would  permit  the  prediction  of  the  fatigue 
life  of  the  associated  panel  structure.  Dowell  further 
suggested  that  conventional  fatigue  data  could  be  used 
for  estimating  the  fatigue  life  of  the  panel  structure 
specifically  at  flutter  conditions.  Xue,  et.al  studied  the 
fatigue  life  of  beams  with  thermal  effects^.  Frequen¬ 
cies  and  stresses  for  the  limit-cycle  behavior  of  these 


The  daniagc-dcpendcnt  stifTness  components  ([A'], 
[iy]i  and  [-0^])  arc  described  in  terms  of  the  damage- 
dependent,  angle-ply  stiffness,  [QT*  Talcrja  defines 
a  parameter,  to  describe  damage  within  \Q'].  The 
stiffness  of  the  composite  laminate  with  damage  is 
found  as 

=  I  Q'ii{l,z,z^)dz  (3.2) 

where 

The  properties  for  the  principal  directions  of  a  lamina 
are  given  by 

Ei  =  E°  +  2^[C^  + 

E2  =  E°  +  2C{Cs  +  -  Ci6«/2°i) 

1^12  =  1^12  +  ^ 

Gi2  =  G°2  ~  (3-3) 

where  j  ^12>  ^12  are  properties  of  the  un¬ 

damaged  lamina.  C3,  Cg,  C13,  and  Cie  are  phenomeno¬ 
logical  constants.  The  damage  parameter,  represents 
the  crack  density,  the  magnitude  of  which  is  subject 
to  a  damage  growth  law.  The  damage  growth  law  for 
matrix  cracking  is  given  as^° 

doc22  =  ^^kG^dN  (3.4) 

where  a22  is  the  ISV  representing  matrix  cracking,  k 
and  n  are  material  constants,  G  is  the  strain  energy 
release  rate,  S  is  the  crack  surface,  and  N  is  the 
number  of  cycles.  The  ISV  is  converted  to  the  damage 
parameter, 

Game  Theoretic  Control 

The  presence  of  damage  to  the  primary  structure,  or 
the  partial  loss  of  distributed  actuator  systems,  presents 
unique  difficulties  in  feedback  control.  Hence,  a  game- 
theoretic  control  methodology  is  derived  in  this  paper 
for  distributed  parameter  systems  subject  to  model  un¬ 
certainties  arising  from  classes  of  damage.  In  some 
recent  work,  conventional  linear  quadratic  regulator 
(LQR)  and  linear  quadratic  Guassian  (LQG)  control 
strategies  have  been  employed  to  enhance  vibration 
suppression  in  aeroelastic  problems^^.  However,  it  is 
well  known  that  LQG  control  designs  can  exhibit  ar¬ 
bitrarily  poor  response  for  unmodelled  perturbations, 
and  adaptive  materials  can  exhibit  different  response 
characteristics  over  time  due  to  evolution  of  their  con¬ 
stitutive  laws.  Hence,  innovative  strategies  for  the  re¬ 
alization  of  controllers  associated  with  game-theoretic 


control,  appropriate  for  direct  treatment  of  the  govern- 
ing  partial  dlfTcrcntial  equations,  arc  derived.  Game- 
theoretic  control  accounts  for  a  class  of  disturbances 
subject  to  well-defined  classes  of  initial  conditions.  Our 
results  show  a  significant  improvement  compared  to  tra¬ 
ditional  LQR/ LQG  methods.  Although  some  efforts 
have  been  made  to  derive  active  control  strategics  for 
aeroelastic  structures  using  smart  materials,  this  re¬ 
search  is  in  the  early  stage  of  development.  The  authors 
are  aware  of  no  research  that  treats  the  fundamental  is¬ 
sue  of  robustness  of  control  design  with  respect  to  actu¬ 
ator  degradation  expected  during  routine  operation,  let 
alone  the  possibility  of  sensor  or  actuator  failure.  While 
numerous  research  papers  on  the  use  of  smart  materi¬ 
als  for  active  control  are  available,  the  works  of  Scott 
and  Rogers  are  particularly  pcrtinent^Liz  Scott,  an 
LQR  strategy  is  employed  to  design  piezoelectric  con¬ 
trols  to  reduce  panel  flutter  and  raise  the  value  of  dy¬ 
namic  pressure  associated  with  flutter.  The  analytical 
results  are  e.xtremely  promising.  On  the  other  hand, 
the  work  of  Rogers  uses  shape  memory  alloys  to  reduce 
crack  growth  in  composites.  One  potential  drawback 
with  the  approach  taken  by  Scott  is  that  the  LQR  de¬ 
sign  strategy  does  not  account  apriori  for  injected  dis¬ 
turbances  or  material  change  which  will  develop  in  fa¬ 
tigue  or  damage.  Likewise,  Rogers  clearly  demonstrates 
that  shape  memory  alloys  have  the  potential  to  sig¬ 
nificantly  reduce  local  stress  concentrations;  however, 
an  active  control  strategy  based  upon  some  optimality 
principle  does  not  exist.  Furthermore,  the  approach  in 
Rogers  deals  with  degradation  of  the  composite  lami¬ 
nate  while  the  integrity  of  the  control  is  unaffected. 

To  motivate  the  design  of  control  strategies  for  the 
aeroelastic  model  with  damage  accumulation,  the  gov¬ 
erning  equations  can  be  interpreted  as  a  two  player 
game.  The  simplest  explanation  is  achieved  by  exam¬ 
ining  the  prior  governing  equations,  but  including  con¬ 
trol.  In  this  case,  the  two  right  hand  side  terms  have 
“adversarial”  roles, 

[[ii:o]  +  [A]]M  +  [B]M  (4.1) 
The  stiffness  damage  and  aeroelastic  terms,  [[iifi)]  + 
[A]],  tend  to  destabilize  the  system,  while  the  con¬ 
trol  terms,  [H],  are  selected  to  offset  this  effect.  The 
destabilizing  terms  represent  localized  damage  from  mi¬ 
crocracking  or  delamination,  actuation  loss,  and  other 
unraodelled  disturbances.  The  control  design  philoso¬ 
phy  employed  in  this  paper  can  be  formulated  by  first 
considering  the  abstract  Cauchy  problem  on  a  Hilbert 
space  H 

x{i)  =  Ax{t)  -f  Bu{i)  (4.2) 

x(0)  =  Xq 

y{x)  =  Cx{i) 

where  B  Q  jEr)(i.e.,  the  linear  operator  from  U 
into  H)j  Cf  is  a  Hilbert  space,  and  A  is  the  infinitesimal 


following  modification^'^  of  a  result  due  to  Zabezyk^^’^^ 
is  required. 

Lemma  4.1  Suji-posc  that  K  G  and  (C,  A)  is 

dcicctahlc.  If  there  is  an  operator  II  >  0  satisfying  the 
fjiapunov  equation 

(IIx,  Ay)  H-  {Ax,  Uy)  —  -(Cx,  Cy)  -  (/v  x,  Ky)  (4.13) 

Vx,  y  G  dornA  then  A  is  stable.  Now,  the  following  re¬ 
sult  can  be  derived,  in  a  manner  similar  to  the  approach 
in  Rhec^^  and  Fabiano^'^. 

Theorem  4.2  Suppose  that  (C,  A)  is  detectable.  If  11 
is  a  nonnegative  definite  solution  of  (4^13),  then  the 
feedback  control  u{i)  ~  ~  D*  Y[x{t)  stabilizes  (4-5) 
for  all\\LA\\  <7,  11^-bII  <  7- 

Proof;  In  Hilbert  space,  it  must  be  shown  that  the 
closed  loop  operator  Ac  ~  A-f  AA  — 
is  stable.  Note  that  (4.12)  can  be  written  as 

(fix,  Ay)-h(Ax,  Hy)  4-  (Qx,  y)~ 

-h  FF^])nx,  Uy)  =  0  (4.14) 

for  all  a:,  y  G  domA.  It  follows  that 

(fix,  Aev)  -f  (AcX,  Hy)  -  (Hx,  Ay)  -f  (Hx,  A  Ay) 

-  {Uxr{B  4-  AB)R-^B^Uy) 

4-  (Ax,  riy)  4-  (A Ax,  Hy)  4-  {Qx,  y) 

-  {{BR-'^B*  -  7^[£)jP*  4-  FF*])nx,  Hy) 

+  FF*)Ux,Uy)  -  {Qx,y) 
-{ABR-^B^Ux,Uy) 

The  underlined  terms  cancel  to  zero,  so  that 
(Hx,  A^y)  +  {AcX,  Ily)  = 

-  (Cx,  Cy)  -  D*ny)  -  7=*(F*nx.  F'Hy) 

-  (B’Hx,  R-^B^Hy)  +  (Hx,  -  {Ex,  Ey) 

+  {DLaEx,  Hy)  -  (Hx,  FLBGir'^B''Uy) 

-  {FLBGR-'-B*nx,'n.y) 

Continuing,  we  have 

(IIx,  Acy)  +  {AcX,  ny}  =  -{Cx,  Cy) 

-T^({P*n  -  ^LaE)x,  (jD*n  -  ^iAE)y) 

-  +  LBGR-^B*)nx, 

{-i^F" +  LBGR-^B*)U.y) 

-  {RR-^B-YLx,  R-'-B*By)  -  {Ex,  {I  -  ^L\LA)Ey) 

-  {GR-'^B'Iix,  {I  -  ^L'BLB)GR-^B'ny) 

Now  if  (IL^II  <  j  and  HIbII  <  7,  where  ll.||  is  the 
operator  norm,  then  the  right  hand  side  can  be  written 


as-(Cx,Cy)-(i^x,/^y)  ,  so  that  the  result  follows  from 
Lemma  3.1.  One  should  carefully  note  tliat  the  above 
proof  relies  on  the  fact  that  DD^:  II  — >  //,  for  example, 
0*  :II  II^.  The  proof  in  Rhce^'*^  and  Fabiano^'^  is 
structured  such  that  O'  :  11  II. 


Control  Formulation 


As  an  illustration,  wc  examine  a  bcam(Fig.l)  prob¬ 
lem  with  piezoelectric  actuation.  In  practice,  the  model 
uncertainty  may  be  due  to  imprecise  actuator  charac¬ 
terization,  actuator  loss  of  authority  due  to  material 
degradation,  mis-calibration,  microcracking  or  dclami- 
nation.  In  the  absence  of  model  uncertainty,  the  gov¬ 
erning  equations  characterizing  the  dynamic  response 
due  to  piezoelectric  actuation  can  be  idealized  as 


-  9^[H{F)\u{i) 


(5.1) 


where 


H{x)  =  H{x  -a)-  II(x  -  /3)  (5,2) 


is  the  Heaviside  function  and  p  is  the  material  density,  h 
is  the  laminae  thickness,  EI{x)  is  the  bending  modulus. 
Cl  is  the  viscous  damping  coefficient,  ki  is  the  stiffness 
constant  of  the  beam  foundation,  and  g  is  the  piezoelec¬ 
tric  gain.  To  account  for  several  regions  of  degradation 
or  damage,  structured  uncertainty  is  represented  in  the 
system  as 

+  (5.3) 

The  boundary  conditions  are  given  by 

!/(i.O)  =  yr(i|0)  =  0 

[(£;i(x)  +  €£;iABl(®))Tii(t,x)]^_,  =  0 

A[(£;j(x)  +  CBxAsx(x))y„(t,x)]^^,  =  0 


and  initial  conditions  are 


y(0,  x)  =  yo(x),  ye(0,  x)  =  yi(x) 

In  this  equation  y(t,  x)  represents  the  transverse  dis¬ 
placement  of  the  beam  at  time  i  and  position  x  along 
the  beam  of  length  L.  This  system  can  be  formulated 


bounds  with  respect  to  the  difTerent  uncertainty.  In  fig¬ 
ure  14,  the  distribution  of  EI{x)  from  the  rnicrocrack 
density  data  derived  in  this  paper  is  used  in  conjunction 
with  the  game  theoretic  control  law.  Figure  15  shows 
the  control  influence  matrix  5(x).  In  this  case,  the  max¬ 
imum  allowable  perturbation  has  been  selected  as  the 
maximum  from  the  rnicrocrjick  density  over  a  portion  of 
the  plate.  From  these  performance  results,  we  can  see 
the  game  theoretic  control  witlrthe  multi-player  game 
theory  is  more  robust  than  the  LQR  control  law. 

CONCLUDING  REMARKS 

The  use  of  composites  and  active  mcitcrials  requires 
the  characterization  of  naturally-occuring,  unavoidable 
damage  effects.  Dynamic  and  aeroelastic  response  is  de¬ 
pendent  upon  damage.  In  turn,  the  growth  and  distri¬ 
bution  of  damage  is  not  uniform;  rather,  damage  growth 
is  dependent  upon  the  presence  of  the  aerodynamic 
loads.  Modal  and  stability  characteristics  provide  a 
nonevasive  measure  for  damage  detection  and  system 
identification.  Active  control  of  aeroelastic  structures 
must  account  for  failure  of  the  primary  structure,  sen¬ 
sors,  and  actuators. 
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Fig.  I .  Aeroclastic  stability  and  control  is  predicted  for  a  composite  wing  structure 
with  distributed  damage.  I  lie  left  view  illustrates  the  wing  structure. 

The  right  view  illustrates  the  "beam"  case  of  panel  flutter. 


Fig.  2.  Damage  is  not  uniform  and  depends  upon  the  aeroelastic  response. 
0  The  distribution  of  damage  is  shown  for  an  outer  ply  of  the  wing  structure 

( left  view  -  early  damage;  right  view  -  prior  to  failure ). 


ig.  3  The  frequency  and  damping  characteristics  Fig.  4.  Aeroelastic  instabilities  change  with  time  (N) 
depend  on  the  dynamic  pressure.  Flutter  occurs  as  two  due  to  the  grorvth  of  damage.  The  flutter  boundary 
(or  more)  modes  coalesce  to  a  common  frequency  and  reflects  the  spatial  and  temporal  distribution  of  stress 
the  system  becomes  undamped.  and  damage. 
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Fig.  11.  Initial  Response  at  Center 
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Fig.  12.  Variation  of  El(x)  (gamma=.0064) 
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Fig.  13  Initial  Response  at  Center 
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Fig.  14  Variation  of  El{x)  (using  crack  density  data) 
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